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ACCURACY OF AERODYNAMIC-HEATING PREDICTIONS 

By A. L. Nagel and R. A. Hanks 
Boeing Airplane Company 


INTRODUCTION 


It has been shown in previous discussions that aerodynamic heating 
has a major role in determining the performance and safety of the Dyna- 
Soar vehicle. A careful examination of the methods which have been 
used to calculate aerodynamic heating rates during reentry, and a com- 
parison of those same methods with test data is a necessary part of 
performance evaluation. 
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sweep angle 
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p. viscosity 

p density 

Subscripts: 

D fraction in dissociation; based on diameter 

EXP experimental 

L 
1 
1 
1 
5 

TH theoretical value 

w evaluated at wall temperature 

oo free-stream value 

o stagnation condition 


e exterior condition 

r recovery 

s stagnation point or line value 


DISCUSSION 


The most severe heating rates on winged hypersonic vehicles will 
occur at the nose and leading edges. The areas involved are relatively 
small, however, and may admit structural solutions (local cooling or 
refractories) which are not practical for the remainder of the vehicle. 
In such a case, the lower- surface material at its most forward point 
may also become a critical heating point, and may be nearer its tempera- 
ture limit than either the nose or the leading edge. Other points which 
would have high local heating rates would be protrusions below the lower 
surface, such as ventral fins or a dihedral ridge line. Some early 
Dyna-Soar configurations had such features and were eliminated for that 
reason. By the end of the Phase I studies, both the Boeing and the 
Martin-Bell teams had arrived at configurations having simple geometry 
in regions of high heat transfer. Upper surface complications are less 
important as the overall heating level for the upper surface is very 
much lower than for the lower surface . The present Dyna-Soar configura- 
tion has only four critical heating points. These are indicated in 
figure 1 as the nose, the leading edge, the lower surface just aft of 
the nose, and the dorsal leading edge of the fin, which is critical at 
low angles of attack. 
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Complete simulation of the reentry environment is not possible in 
any of the ground facilities which must be used to provide the bulk of 
heat-transfer test data. For this reason extrapolation of the test 
results to the flight condition by some theoretical method is necessary. 
The study of prediction accuracy cannot be limited to an examination of* 
the data scatter, but must include an evaluation of the theoretical 
method as well. The combination of a rigorous theoretical approach and 
test data taken in facilities which simulate the important aspects of 
the flight environment allows a high degree of confidence in the pre- 
diction. In other cases, the theory may be too idealized to lend cre- 
L dence to extrapolations. 

1 
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1 Heat-Transfer Equations 

5 

The equations used for calculating laminar heating rates both for 

he reentry condition and for the following comparisons with test data 
are : 
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The constant in equation (l) is 0.793 for the axisymmetric stagnation 
point and 0.576 for the two-dimensional stagnation line. Equation (2) 
is used for calculating heat-transfer rates away from the stagnation 
point. Both equations (l) and (2) are from the work of Kemp, Rose, 
and Detra, (ref. l) which is an extension of the earlier work of Fay 
and Riddell (ref. 2). These equations were selected as a basis for 
reentry-heating calculations because they are the most rigorously 
developed methods available, and because they are in good agreement 
with the test data, as will be shown. The expressions were originally 
obtained by numerically integrating the boundary- layer equations, using 
the real-gas equation of state and the Sutherland viscosity law/ The 
cases specifically considered corresponded to the axisymmetric stagna- 
tion point, the unswept stagnation line, the flat plate, and a limiting 
pressure gradient case. Applying simple geometric corrections for flow 
pattern allows the results to be used for swept leading edges as well. 
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Calculations by Beckwith (ref. 3) have further shown that equa- 
tion (l) also results from ideal-gas calculations for both the swept 
and unswept leading edges except that the Lewis number term is, of 
course, missing. 


The velocity gradient used in applying equation (l) was based on 
the modified Newtonian pressure distribution, which is within a few 
percent of the best known values. Equilibrium dissociation was assumed 
and the values of viscosity calculated by Hansen (ref. 4) were used, and 
the Lewis number was taken as 1.4. Use of the higher viscosity values 
has been found to improve agreement with test data. Evaluated in this ■ L 

way heat-transfer rates are 5 4o 45 percent higher than those obtained 1 

by the method of Fay and Riddell in reference 2. 1 

1 

Since the calculations on which equation (l) is based assumed the 3 


Sutherland law for viscosity, use of another viscosity law might appear 
to invalidate the equation. Recent unpublished calculations of Beckwith 
and Cohen at the Langley Research Center have shown, however, that the 
form of the equation does not depend on the viscosity law or even upon 
the equation of state. It appears then, that use of equation (l) with 
the best available fluid properties will provide the best estimate of 
heat transfer. 

In the form shown here, the heat-transfer distribution function 
(eq. (2)) depends only on the local pressure and flow velocity. Simpli- 
fying assumptions first suggested by L. Lqes (ref. 5) are required to 
eliminate the dependence on the local transport properties and the 
pressure gradient. The simplification has been found to be satisfac- 
tory for shapes without sharp corners, such as the present Dyna-Soar 
nose . 

In the original development, equation (2) was intended for appli- 
cation to two-dimensional or axisymmetric bodies at zero angle of 
attack. Application to less simple shapes can be accomplished by 
replacing the radius terms with an equivalent radius which expresses 
both the body shape and the streamline pattern which occurs on it. 

Equations (l) and (2) supply the required laminar flow heating- 
rate estimates at all the critical points. At the nose and forward 
lower surface the Reynolds number is in the range for which laminar- 
flow heating rates are higher than the turbulent flow rates. The 
leading edge, however, is limiting in turbulent flow at velocities 
less than about 19,000 ft/sec. It might appear at first that turbu- 
lent flow cannot exist at the leading-edge stagnation line. This is 
not true for the swept leading edge, as the flow along the stagnation 
line can become turbulent. The possibility of turbulent low leading- 
edge heating rates must therefore be considered in aerodynamic heating 
calculations. The tendency of the boundary- layer secondary flow to 



promote transition at very low Reynolds numbers makes it especially 
important to consider turbulent boundary- layer flow. The expression 
used to calculate turbulent leading-edge heat transfer for the reentry 
condition as well as for the comparisons with wind-tunnel data which 
follow is: 



This expression was developed by Beckwith and Gallagher (ref. 6). A 
similar expression can be obtained by applying geometric corrections 
to turbulent-f low flat-plate theory. 


Experimental Comparisons 

Stagnation point .- The experimental data for hemispherical stagna- 
tion point available for comparison with the theory are presented in 
figure 2. Data are taken from shock- tube experiments (ref. 7), wind- 
tunnel tests (ref. 8), and free-f light tests (ref. 9) and have been 
divided by the theoretical value for the same conditions. A ratio 
of 1.00 therefore indicates perfect correlation of theory and experi- 
ment. The data are shown to scatter from O.65 to 1.4 times the theory 
with the average very nearly 1.0 over the entire velocity range. It 
is believed that this large scatter reflects experimental errors rather 
than fluctuations in the actual heat-transfer rates. This view is sup- 
ported by the random nature of the scatter. 

The comparison of figure 3 lends further support to this explana- 
tion of the scatter. Heat-transfer data from an Atlas (ref. 10) reentry 
flight are compared with the theory as a function of time. The theory 
and experimental curves are of similar shape and with almost identical 
peak values, but with an offset of about 2 seconds in time. If plotted 
in figure 2, these data would have shown a scatter of about 50 percent 
around the theoretical curve. The heat-transfer rates shown in figure 3 
are calculated from the temperature response of the skin as recorded by 
thermocouples installed in plugs in the skin. The same characteristic 
lag of experiment behind theory was observed in many flights. After 
several other explanations had been ruled out by the consistency of the 
lag, similar thermocouple installations were calibrated in ground tests. 
Lags in heating rate were found to occur which, when extrapolated to the 
flight conditions, are of the same order as those observed in flight. 

£hey cannot be said to be precisely the same, as no two thermocouple 
installations showed exactly the same lag. However, the thermocouple 
lag does seem to provide a reasonable explanation of the offset. Another 
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point worth mentioning is in the comparison of the theory of reference 2 
with these same data, which is shown to "be approximately 10 percent too 
low at peak heating. Equation (l) was originally selected in preference 
to the Fay and Riddell theory on the basis of the shock-tube data pre- 
sented in figure 2. Although both theories fell well within the scatter 
of the data, the values from the Fay and Riddell method were about 
10 percent below the average of the data at high speeds. Re-evaluating 
the theory to improve the agreement with the average of the shock-tube 
data brings it into almost exact agreement with the flight data. 


Hot-gas radiation heating and variation in wall catalytic effect L 
are two effects not reflected in these data. Calculations based on best 1 
available information show that the radiative heat -transfer rate is very 1 
small compared to the convective rate for nose radii of 1 foot or less. 1 
This result is consistent with the conclusions of previous papers in 5 
this conference. The effect of wall catalytic effect may be signifi- 


cantly favorable if nose coating materials can be developed which do not 
catalyze recombination at the wall. Reductions of over 50 percent are 
theoretically possible. Neglecting this effect is conservative, and 
appears to be most realistic at present, as noncatalytic materials have 
not been developed. 

Laminar distribution .- Nose hemisphere and afterbody data are com- 
pared with the theory in figure 4. The agreement is shown to be satis- 
factory for cones (ref. ll) . Recent Langley data on delta wings at 
angle of attack presented by Bertram and others at this conference are 
also in good agreement. In making these comparisons equation (2) has 
been used with an effective radius to correct for nonaxi symmetric 
shape in the manner previously suggested. This correction requires a 
knowledge of the streamline patterns as well as the shape of the body. 

For the cone, the streamline pattern is based on values from the Kopal 
tables (refs. 12 and 13). For the delta wing, streamline patterns were 
calculated from a correlation previously made of oil flow patterns 
obtained in delta wing tests. 

Laminar leading edge .- Laminar leading edge data from several wind- 
tunnel tests (refs. 14 to l6 and unpublished Boeing and Langley test 
data) are presented in figure 5. The agreement is shown to be excellent 
over the entire range of sweep angles. None of the experimental leading- 
edge data shown are in the total temperature range for which real gas 
effects would be distinguishable. However, the theoretical expression 
indicates that for the highly swept leading edge the real-gas effect is 
much smaller than at the stagnation point, and even for that case the 
predicted (and experimentally supported) effect is only 15 percent at 
the velocity for maximum reentry heating. 

Turbulent leading edge .- Experimental data (ref. 6 and unpublished 
Boeing test data) for turbulent stagnation-line heat transfer are com- 
pared to theoretical values in figure 6. The data are predominantly 
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below the theory, indicating principally the difficulty of obtaining 
turbulent flow on the leading edge. The low points at 10°, 20°, and 75° 
sweep angles are apparently transitional. The low data points taken 
on delta wing at angles of attack of 30° and 3k° are affected by the 
presence of the wing, which distorts the inviscid flow field. The rest 
of the data are in very good agreement with the theory. 

Extrapolation to the flight condition is still somewhat uncertain 
because the theory has assumed ideal gas relations throughout. Some 
information regarding the validity of ideal gas heat-transfer calcula- 
tions in a real gas environment is afforded by an examination of the 
effect in laminar flow, for which the theory is well developed. Com- 
parisons of ideal-gas solutions for both the stagnation point and the zero 
pressure gradient flat plate have been found to be within 10 percent of 
rigorous real-gas solutions at speeds up to satellite velocity. A more 
significant comparison is presented in figure 7. Experimental turbulent 
heat— transfer data (refs. 10 and 17 to 20) in the real— gas temperature 
range are compared with the ideal-gas reference temperature method. The 
agreement of the theory and experiment is very good over the entire veloc- 
ity range. The good agreement between the normalized heat transfer in 
free flight and in the shock tube tends to eliminate the possibility of 
fortuitous agreement. This agreement also indicates that the shock tube 
can provide fundamental information about real-gas effects in turbulent 
flow just as it has been used in the past to study laminar flow stagnation 
point heat transfer. 


CONCLUDING REMARKS 


It has been shown that there exists a unified theoretical method 
for laminar-flow heat transfer which is applicable to critical tempera- 
ture locations on the Ifyna-Soar vehicle. The method rigorously includes 
real-gas behavior and other phenomena having significant effects on 
heat transfer . Test data have been presented which confirm the ability 
of the method to account for such effects over a wide range of conditions. 

In the turbulent— flow case no similar well developed theory exists. 
There is, however, a compensation for this deficiency in the relative 
insensitivity of the turbulent boundary layer to any influence other 
than local pressure and velocity. Some turbulent -flow heat- transfer 
data in the speed range corresponding to reentry maximum heating do 
exist, and these data are in agreement with semi empirical methods now 
in use . 

From these comparisons it appears that existing methods will sat- 
isfactorily predict aerodynamic heating during reentry for the critical 
locations on the present configuration. Further testing is required to 
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substantiate this conclusion with the emphasis on data for the con- 
figuration specifically chosen. Further testing is also desirable to 
reduce the uncertainty caused by scatter in available data. These data 
appear to reflect experimental errors, rather than fluctuations in 
actual heating rates, so that a design based on these data alone would 
incorporate unneccessarily large margins in temperature capability, 
with corresponding weight and performance penalties. 

The methods used for theoretical calculations can be extended to 
other locations on the vehicle. As previously stated, the accuracy of 
calculations over the rest of the vehicle is less important, as one of 
the four points discussed will always be nearer its l i . rn . i t temperature. 
Future alterations of configuration or materials may cause other points 
to become critical. 
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DYNA-SOAR REENTRY VEHICLE 

CRITICAL HEATING LOCATIONS 
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Figure 2 



STAGNATION-POINT HEAT TRANSFER 
FLIGHT TEST DATA 
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NOSE -AFTERBODY HEAT TRANSFER 
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LEADING-EDGE LAMINAR FLOW 
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Figure 5 


LEADING-EDGE TURBULENT FLOW 

O SWEPT CYLINDER; R^ D » I0 6 ; REF. 6 

WING L.E.; A = 75°; R^ D = 0.75 X I0 6 ; UNPUBLISHED 
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Figure 6 


TURBULENT HEAT TRANSFER 
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Figure 7 


